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We present the first model independent search for three-jet hadronic resonances within multijet 
events in ^/s — 1.96 TeV pp collisions at the Fermilab Tevatron using the CDF II detector. Pair 
production of supersymmetric gluinos and squarks with hadronic R-parity violating decays is em- 
ployed as an example of a new physics benchmark for this signature. Selection criteria based on 
the kinematic properties of an ensemble of jet combinations within each event help to extract sig- 
nal from copious QCD background. Our background estimates include all-hadronic tt decays that 
have a signature similar to the signal. No significant excess outside the top quark mass window is 
observed in data with an integrated luminosity of 3.2 fb~^. We place 95% confidence level limits on 
the production cross section i7{pp — ^ XX ) x BR(gr(/— > 3 jet -I- 3 jet) where X, X' = g, q, or q, with 
q, q ^ g -\- jet, as a function of gluino mass, in the range of 77 GeV/c^ to 240 GeV/c'^. 

PACS numbers: 13.85.-t, ll.30.Pb 



Most searches for new physics at high energy hadron 
colliders use signatures that require leptons, photons, or 
missing transverse energy (.^t) [H in order to suppress 
backgrounds from QCD. Final states with multijets and 
.^T have also been explored 0, Q • 

In this letter, we present a first new physics search in 
an entirely hadronic channel with no .^t signature us- 
ing data collected with the Collider Detector at Fermilab 
(CDF). This data set corresponds to an integrated lumi- 
nosity of 3.2 fb~^ of pp collisions at ^/s = 1.96 TeV at 
the Tevatron collider. The search utilizes a novel ap- 
proach [1, 1^: an ensemble of all possible jet triplets 
within an event consisting of at least six jets is used 



to extract a signal from the multijet QCD backgrounds. 
We model the possible new physics origin for this signa- 
ture with pair production of SU{i)c adjoint Majorana 
fermions each one decaying into three quarks 0| ■ This 
search is sensitive to models such as hadronic R-parity 
violating supersymmetry (RPV SUSY) [^] with a gluino, 
chargino, or neutralino lightest superpartner, as well as 
the hadronic decay modes of pairs of top quarks or fourth 
generation quarks and it complements existing di-jet res- 
onances searches at hadron colliders. Moreover, it does 
not require any 6-quark jet identification which is an im- 
portant tool, often used for top quark identification. 

The CDF II detector is a multi-purpose particle de- 
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tector consisting of tracking and calorimeter systems [8| . 
The data were collected using an online event selection 
that requires at least four calorimeter jets Q with uncor- 
rected transverse energy Et > 15 GeV. A jet is formed 
by a cluster of calorimeter towers and reconstructed with 
a cone algorithm using a fixed cone of AR = 0.4 [lo| . 
with AR = y/Ar]^ + ^ 

In the online selection an 
additional request is made for the sum of the transverse 
energy of all clusters to be larger than 175 GeV. At the 
analysis level, jet energies are corrected to account for ef- 
fects such as non-linearities in the_detector response and 
multiple pp collisions in an event 



PYTHIA RPVGIuino 



11|. 



Events are selected with at least six jets with transverse 
momentum (px) greater than 15 GeV/c and \r]\ < 2.5. 
The scalar sum of the most energetic six jets' pt, '^Pt, 
is required to be greater than 250 GeV/c and events with 
> 50 GeV are removed. Multiple interactions, result- 
ing in the reconstruction of more than one primary vertex 
in the same event, contribute to the multijet background. 
We require at least one primary vertex and discard events 
with more than four primary vertices. To further reduce 
this background, we require jets in an event to originate 
from near the same point on the beamline. We asso- 
ciate tracks with each jet where possible [l^l by requir- 
ing AR between the track and the jet to be less than 
0.4. The mean z-coordinate of all tracks associated with 
each jet {zj for the j^^ jet), and the associated stan- 
dard deviation {5{zj)) are determined. Events with jets 



that have 



> 60 cm are discarded. We then evaluate 



the standard deviation of the Zj of all jets in the event 
((5(zg^jj)) and select events that have at least four jets 
with 5{zj) < 4 cm, and 5{zg^^^ < 0.5 cm, consistent with 
the resolution of tracks associated with jets. Once the se- 
lection is applied, pileup effects are significantly reduced. 
Since we select events with at least six jets, we consider 
an ensemble of 20 (or more) possible jet triplets. We dis- 
card those triplets that have more than one jet with no z 
information. In addition, all jets in the triplet must have 
5{zj) < 2.5 cm, and originate from within 10 cm of the 
primary vertex of the event. 

The biggest challenge of this analysis is to reduce the 
large multijet QCD background. To extract signal from 
this background, we apply the following technique: for 
every accepted triplet we calculate the invariant mass. 

Triplets made of un- 



PT, while sig- 



Mjjj, and scalar sum px, J^jjjPt 
correlated jets tend to have MjjjC 

nal triplets should have Mjjj as close to the mass of the 
decaying particle as allowed by jet energy resolution. We 
then select triplets with Pt — ^jjjC > A, A being a 
diagonal offset as illustrated in Fig.[TJ The diagonal offset 
values are optimized for the best signal over background 
ratio separately for each hadronic resonance mass in this 
search. The optimized diagonal offset selection greatly 
reduces the QCD background and the contribution from 
incorrect combinations of jets. We note that for a small 
fraction of events it is possible for multiple triplets to 
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FIG. 1. Distribution of Mjjj versus 'YlijjjPT for a pair- 
produced RPV gluino with invariant mass 190 GeV /c^ gen- 
erated with PYTHIA MC. Triplets to the right of a diagonal 
offset C^ZjjjPT — MjjjC = A), indicated by the dashed line, 
are are kept. The inset shows the Mjjj distribution for the 
RPV signal MC and with no QCD background after a diago- 
nal offset of 195 GeV/c along with a Gaussian plus a Landau 
fit; the Landau shows the combinatorial contribution within 
the signal jet ensemble. The QCD background distribution re- 
sembles that of the combinatorial contribution, because they 
are both due to effectively uncorrelated triplets. 



pass all selection criteria. 

The QCD background is estimated from a 5-jet data 
sample, which is statistically independent of the signal 
sample of > 6 jets (for brevity referred to as 6-jet). The 
5-jet Mjjj distribution is rescaled by the ratio of the 6- 
jet to 5-jet population in each X^jjj Pt bin. A Landau 
function is chosen [3] to fit the scaled 5-jct Mjjj distribu- 
tion. The Landau parameters extracted from the scaled 
5-jet Mjjj distribution vary by less than 2 GeV/c^ from 
similar fits to the 6-jet sample, indicating that the scaled 
5-jet sample describes the background in the 6-jet sam- 
ple well. The contribution to the background from tt pair 
production is estimated using the pythia Monte Carlo 
(MC) gejnerator [l3| followed by the CDF detector simu- 
lation !14|. These events were generated assuming a top 
quark mass of 172.5 GeV/c^ and production cross section 
of 7.5 pb. To ensure a proper fit to the QCD background, 
the fit is blinded to the mass region corresponding to the 
top quark, 153 GeV/c2< Mjjj <189 GeV/c^. Addition- 
ally, we find that truncating the Landau fit for lower val- 
ues of A gives an improved description of the QCD back- 
ground. The Landau parameters extracted from the fits 
vary smoothly as functions of the diagonal offset value. 
We now have a firm prediction for the QCD background 
and fix the parameters when we fit for signal. 

The signal is modeled using the pythia MC generator. 
The process pp — > XX where X, X' = g, q, or q is 
simulated at several gluino mass values, ranging from 
74 GeV/c2 to 245 GeV/c^ with hadronic uds RPV SUSY 
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decays turned on, allowing gluino decays to three light 
jets. Two scenarios of squark masses are considered (0.5 



TeV/c2 < m~< 0.7 TeV/c^ 



10 GeV/c2) and o 



9 ' q g 

were found to give equivalent acceptances. 

The acceptance of the trigger, reconstruction, and se- 
lection requirements for signal events is determined by 
fitting the pair produced RPV gluino MC with a Lan- 
dau plus Gaussian function, corresponding to the com- 
binatorial contribution and signal peak respectively. An 
example is shown in the inset of Fig. [TJ The Gaussian is 
integrated in a ztltr range to extract the number of signal 
triplets. This procedure is repeated for various diagonal 
offset values and the optimal offset for each hadronic res- 
onance mass is determined. The acceptance, calculated 
for these optimal offset values, is 5x 10~^, constant within 
20% across all gluino mass points. 

The expected sensitivity of this analysis in the absence 
of signal is determined with a set of background-only 
experiments (pseudoexperiments) . A pseudoexperiment 
is constructed with the background modeled by a Lan- 
dau function whose parameters are chosen randomly from 
within the range allowed by the background shape fits, 
with the expected amount of tt added. Each pseudoex- 
periment is fit with the Landau background shape pa- 
rameters fixed, and a signal Gaussian whose position is 
determined by the mass point being fit, and whose am- 
plitude and width are allowed to vary within a range 
determined by the expected signal shape. The number 
of signal triplets allowed by each pseudoexperiment is ex- 
tracted by integrating the Gaussian in the same way as 
in the acceptance calculation. 

Two broad categories of systematic uncertainties are 
accounted for in extracting a cross section: uncertain- 
ties in the shape of the Mjjj distribution and uncertain- 
ties in the acceptance of the signal. Shape uncertainties, 
determined from background and signal fits, are incor- 
porated in the pseudoexperiments themselves. Accep- 
tance uncertainties arise from modeling the signal Monte 
Carlo and include effects of initial and final state radi- 
ation |l5| (20%), parton distribution functions (PDFs) 



from CTEOjl6| (10%), jet energy scale [ll| (31%) and 
luminosity [l7| (6%) uncertainties. The overall accep- 
tance uncertainty due to these sources is 38%. 

We search for a hadronic resonance in the data for an 
invariant mass (m) 77 - 240 GeV/c^ in 9 GeV/c^ steps, 
consistent with jet energy resolution. For each mass, jet 
triplets are selected by the optimal diagonal offset value. 
The data Mjjj distribution is fit in exactly the same way 
as the pseudoexperiments. Figure [2] shows the Mjjj dis- 
tribution for m = 112 GeV/c2 and 175 GeV/c^. The 
latter fit shows a noticeable excess consistent in mass 
with the hadronic decay of the top quark. The Gaus- 
sian component of the fit integrated from 165 GeV/c^ 
to 185 GeV/c^, corresponding to a ztlcr window around 
the Gaussian peak, gives 11 ± 5 triplets. The number of 
expected QCD background triplets in the same mass win- 
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FIG. 2. Mjjj distributions in 3.2 fb~^ data fitted to a Landau 
(parameters are extracted from fits to the scaled 5-jet Mjjj 
distribution) plus a Gaussian at (a) 112 GeV/c^ (optimal di- 
agonal offset value 155 GeV/c) and (b) 175 GeV/c^ (optimal 
diagonal offset value of 190 GeV/c). The fit function in panel 
(b) includes a Gaussian fixed at m — 175 GeV/c^. 



dow from the Landau function is 8 ± 1. The ti contribu- 
tion to background is evaluated using pythia. It is cross- 
checked with higher order tt MC generators ALPGEN (isj 
and MC@NLO .19'], samples that varied the amount of 
initial and final state radiation, as well as samples that 
varied the PDFs within their uncertainties. These stud- 
ies lead us to expect between 0.5 and 1.1 triplets from tt 
production in the aforementioned mass range. We note 
that ~10% of the triplets in the top mass window origi- 
nate from two or more combinations in a jet ensemble of 
a given event, consistent with the pythia tt simulation. 
We evaluate the significance of the excess using the pseu- 
doexperiment method described above, which includes 
systematic uncertainties on signal acceptance as well as 
the shape of the Mjjj distribution. The observed excess 
is 2 standard deviations (2(t) above the prediction. Ad- 
ditional cross-checks, such as requiring one of the jets to 
have originated from a 6-quark, suggest that the excess 
is consistent with coming from top quarks. 

We do not observe a significant deviation from stan- 
dard model backgrounds anywhere in the data. A 
Bayesian approach is used to place 95% confidence level 
limits on a{pp — )■ XX ) x BK(gg^ 3 jet -I- 3 jet) where 
AT, X' — g, q, or g, with q ^ g + jet, versus gluino 
mass, shown in Fig. [31 The largest excess observed is the 
one previously noted located near the top quark mass. 
We find that our background estimate has a 2.3% prob- 
ability of producing such a deviation. Comparisons to 
the theoretical cross section for a{pp — ^ XX ) x BR(gg^ 
3 jet -I- 3 jet) from pythia corrected by a next-to- leading- 
order (NLO) fc-factor calculated using prospino [1^ are 
shown in the dashed and dash-dot lines for two different 
squark mass scenarios. For a decoupled squark mass (0.5 
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FIG. 3. The observed (points) and expected (solid black line) 
95% confidence level limits on the production cross section 
a{pp XX') X BR(gg -5> 3 jet + 3 jet) where X, X' = p, q, 
or q, including systematic uncertainties. The shaded bands 
represent the total uncertainty on the limit. Also shown is 
the model cross section from pythia corrected by an NLO 
fc-factor (dash-dot line for 0.5 TeV/c^ < m~ < 0.7 TeV/c^ 

dashed line for m~— m~+ 10 GeV/c^). 

TeV/c^ < m~ < 0.7 TeV/c^) we exclude gluinos below 
a mass of 144 GcV/c^ (dashed line). In the case of a 
squark mass which is nearly degenerate with the gluino 
mass {m~= ™~+ 10 GeV/c^) we exclude gluinos below 
155 GeV/c^ (dash-dot hne). 

We have performed a first search for three-jet hadronic 
resonances in a six or more jet final state using a data 
sample with an integrated luminosity of 3.2 fb~^ collected 
by the CDF 11 detector. A novel technique is introduced 
that exploits kinematic features within an ensemble of 
jet combinations that allows us to extract signal from 
the QCD background. We observe no significant excess 
in the data in an invariant mass range from 77 GeV/c^ 
to 240 GeV/c^ and place 95% confidence level limits on 
the production cross section cr{pp — > XX ) x BR(g5— > 
3 jet -I- 3 jet) where AT, X' — g, q, or q, with q, q ^ g + 
jet, versus gluino mass. The results are presented as lim- 
its on RPV gluinos decaying to three jets, but are more 
widely applicable to any new particle with a three-jet de- 
cay mode. Two different squark mass scenarios have been 
considered: decoupled squarks and squarks nearly degen- 
erate in mass with the gluino. We can exclude gluinos 
below 144 GeV/c^ and 155 GeV/c^ respectively. 
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